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The specific kind of diffusion stimulated (rather than suppressed) by the external magnetic field, 
which was predicted for the first time by Schmelcher and Cederbaum in 1992, is considered here for 
the case of high-angular- momentum (i.e., approximately "circular") Rydberg atoms. The coefficient 
of such diffusion is calculated by a purely analytical approach. 



PACS numbers: 05.40.Fb, 05.40.Jc, 51.20.+d, 32.80.Ee 

A well-known feature of the diffusion processes in vari- 
ous gaseous systems (ranging from the laboratory devices 
for plasma confinement to the large-scale astrophysical 
objects) is that it is strongly suppressed by imposition of 
the external magnetic fields. The aim of the present arti- 
cle is to discuss the case when Rydberg atoms, along with 
their other extraordinary properties [e.g., review [lj]), 
may exhibit also a quite unexpected kind of diffusion, 
which is stimulated rather than suppressed by the mag- 
netic field. 

In fact, the possibility of such phenomenon was men- 
tioned for the first time by Schmelcher and Cederbaum as 
early as 1992 @, These authors, using a classical ap- 
proximation for the description of hydrogen-like atom in 
a magnetic field, studied dynamics of a "cyclic" canonical 
coordinate (i.e., the one not entering explicitly into the 
Hamiltonian), such as the center of mass of the atom. As 
a result, they revealed a possibility of chaos development 
in this coordinate and found that, in some circumstances, 
such chaotic motion follows the classical diffusion law. 

Later, a profound mathematical treatment of the same 
problem was given in paper Rs authors used the 
method of separation between the "slow" and "fast" de- 
grees of freedom of a dynamical system. Then, the fast 
chaotic motions were treated as a stochastic force, result- 
ing in the diffusion-like behavior of the slow degrees of 
freedom. 

Experimental observation of such diffusion requires ei- 
ther an extremely strong magnetic field or very weak 
interatomic electric field (i.e., a low binding energy). 
Therefore, it was suggested to observe this phenomenon 
either in compact astrophysical objects with ultrastrong 
magnetic fields (e.g., pulsars) or in the experiments with 
highly-excited (Rydberg) atoms. Unfortunately, none 
of these opportunities was realized by now because, on 
the one hand, the diagnostic possibilities in astrophysi- 
cal studies are too limited and, on the other hand, the 
standard atomic-beam experiments with Rydberg atoms 
are not well suited for studying the diffusion processes. 

Fortunately, the situation changed in the recent 
decade: magneto-optical traps (MOT) became a new 
source of the Rydberg atoms, which are formed due 
to recombination of laser-produced ultracold plasmas 
(e.g., reviews [3, 0]); and such devices are much better 



suited for tracing the long-term diffusional effects. Be- 
sides, quite similar apparatus is used now to create and 
study the antihydrogen by recombining antiprotons and 
positrons [7, 8]. Such antihydrogen atoms are also born 
in the highly-excited states; and strong magnetic fields 
are applied to prevent antiparticles from the escape and 
annihilation. 

However, it should be kept in mind that the above- 
cited theoretical studies were concentrated mostly on the 
case of atoms with zero or very low angular momenta. 
This was quite natural for the old atomic-beam experi- 
ments with Rydberg atoms produced by laser irradiation: 
such atoms can have a huge principal quantum number, 
but their orbital number is always limited to a few units, 
because any absorbed photon can change it only by unity. 
On the other hand, the Rydberg atoms formed by the 
recombination in MOT experiments, in general, should 
have large angular momenta (so that the respective or- 
bital quantum numbers are on the order of the princi- 
pal number). It is not clear in advance if such atoms 
will behave similarly in the external magnetic fields. So, 
a special treatment of the high-orbital-momentum case 
should be performed. Besides, to estimate importance of 
the magnetically-stimulated diffusion in various experi- 
mental setups, it would be desirable to have an analyt- 
ical expression for the diffusion coefficient, with explicit 
dependences on all the relevant parameters. It is just the 
aim of the present work to provide such formula [t| . 

Before proceeding to the detailed quantitative treat- 
ment, let us try to explain qualitatively why the magnetic 
field can stimulate rather than suppress a diffusion. In 
the case of high-angular-momentum (i.e., approximately 
"circular") Rydberg atoms, this can be presented pictori- 
ally in Fig. [TJ Really, if the magnetic field is absent, the 
electron orbit will be exactly closed, and average force 
experienced by the central ion after each revolution of 
the electron will equal exactly zero. On the other hand, 
when an external magnetic field is imposed on the sys- 
tem, then the curvature radius of the electron orbit will 
either increase or decrease under the action of the addi- 
tional Lorentz force. As a result, the orbit will no longer 
be closed; and a position of the electron after each rev- 
olution will be slightly shifted with respect to the initial 
point, as illustrated by the shaded circle in Fig.[TJ There- 
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FIG. 1: Sketch of development of the diffusion-like behavior 
of ion (i) under the action of revolutions of the electron (e) 
perturbed by the external magnetic field. 



fore, the average Coulomb force acting on the ion after 
such revolution becomes nonzero; and this ion will ex- 
perience a kick in a quasi-random direction. A series of 
such kicks is equivalent to a random force, leading to 
the diffusion-like Brownian motion of the ion and, conse- 
quently, of the entire atom. 

Quantitative description of the two-body system, cou- 
pled by Coulomb forces and imbedded in the external 
magnetic field, can be evidently given by the set of equa- 
tions: 

j2 



e dF r e -r, a c dt 



d 2 
dt 2 



.2 r e 



e dr * y R 



(1) 



where m e and rrii are the electron and ion masses, r e 
and vi are their radius vectors, e is the absolute value of 
the electron charge, c is the speed of light, and B is the 
magnetic field vector. 

Next, let us introduce the standard definitions for the 
center of mass of the system 

m e r e + to.^ 



R = 



(2) 



m e + m» 

the relative position of an electron with respect to ion 

r = r e -ri, (3) 

and the following mass parameters: the total mass of the 
system 



M = rrii + m e 



the reduced mass 



rrii m e 



m e 



rrii + m e 
and the relative mass difference 

mi — fTT-e 

v = 



m; 



(4) 
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As a result, the initial set of equations (p} will take the 
form: 
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(R+it)xB. (7b) 



Next, it is convenient to normalize all spatial quantities 
to the initial radius a of the unperturbed orbit: 

r = of , R = oR ; (8) 

and time, to the Keplerian period of revolution: 

t = Tt, where r = 2ny / na 3 /e 2 . (9) 



Then, equations (J7aJ), (|7b|) can be rewritten in terms 
of the dimensionless variables (marked by tildes) as 
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where 



-^f = -(2tt)^ -27r-^^(R+2/f)xb, (10b) 
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Mc rriiC 
is the gyrofrequency of the center of mass, 
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lie 6 m e c 
is the gyrofrequency of the relative motion, 



ui e = 2tt/t = \J e 2 / jia 3 w ^/ e 2 /m e a 3 



(11) 



(12) 



(13) 



is the Keplerian frequency, and b = B/i3 is the unit 
vector in the direction of the magnetic field. 

Since m e <C m„ it can be reasonably assumed that the 
center-of-mass motion is much slower than the relative 
motion of the electron and ion [lfj : 



\dK/dt\ < \dr/dt\ 



(14) 



Consequently, the term dR/dt in the right-hand side of 
Eq. (fl0b|) can be neglected, and this equation becomes 
completely independent of Eq. (|10a[) : 
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Therefore, we can solve this equation alone and then sub- 
stitute the obtained solution to Eq. (jlOal) . 

Rewriting Eq. ((T5|) in the polar coordinate system 
(f, <p) leads to the following set of equations: 
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where dot denotes a derivative with respect to the di- 
mensionless time i. 

Let us seek for the solution of equations (fT!)]) as per- 
turbation of the purely circular motion: 




f + Sf(t) , 

2-Kt + ip + 5<p[t) 



(17) 



where fg = const, tpo = const. 

Then, it can be easily shown that, due to the influ- 
ence of a magnetic field, the mean orbital radius changes 
from 1 (in dimensionless units) to 
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(18) 



while perturbations <5f (t) and 6ip(i) experience the small 
harmonic oscillations with frequency 
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at O e /w e < 1 . (19) 



Not going into details of the perturbed relative motion, 
let us mention only one property, which will be very im- 
portant further: As follows from Eq. (fl~9"|). the electron 
after each revolution will be shifted by some distance ap- 
proximately in the same direction in the course of 
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steps; and then this direction will gradually change, as 
shown by the dotted arrow within the shaded circle in 
Fig. Q] [11]. In other words, n gives the characteristic 
correlation scale of the random walk. 

Now, when the most important features of the effective 
stochastic force in the right-hand side of Eq. (|10a|) are 
established, we can proceed to studying behavior of the 
center of mass R. In general, Eq. (|10al) represents a 
stochastic differential equation, which can be solved by 
various methods. In the simplest approximation, we can 
use the well-known formula from the theory of random 
walks: 



(R 2 (t)) ps AR y/Nt, 



(21) 



Here, AR is the typical length of one uncorrelated step, 
comprising n correlated substeps given by Eq. (|20[) : and 
iVtot is the total number of uncorrelated steps during a 
given time. 

Next, estimating the differential equation (|10a|) at the 
scale of one uncorrelated step, we can easily get the fol- 
lowing estimate: 



AR 



(2tt) — A? At. 



(22) 



where Af and At are the radial and temporal increments 
corresponding to such an uncorrelated step. 

The last-mentioned time interval (in dimensionless 
units) is just the typical number of revolutions in which 
the correlations still survive, which is given by Eq. (|20p : 
so that 
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On the other hand, the above-mentioned radial inter- 
val Af (involving n revolutions of the electron) is just the 
typical size of the shaded circle in Fig [T] Therefore, it 
can be estimated as difference between the average radius 
of the electron orbit in presence of the magnetic field, as 
given by Eq. (TT8]) . and the radius without the field (which 
is equal just to unity). As a result, we get 
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At last, the total number of random-walk steps, ap- 
pearing in Eq. (|2ip , can be evidently obtained as ratio of 
the total time interval to the duration of one step: 



2n e - 

t/n sa t . 



(25) 



Finally, substituting formulas (|22 j) - ([2l)]l into (|2l]) . we 
arrive at 
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which is just the diffusion law, where the coefficient of 
diffusion (in dimensionless units) is 



D 



(27) 



The exact numerical factor appearing here can be hardly 
reliable in the framework of our approximate treatment; 
so we would prefer to drop it out in the final result. 
At last, returning to the dimensional variables, we get 



D 



(28) 



This seems to be the first analytical formula derived for 
the coefficient of magnetically-stimulated diffusion of the 
high-angular-momentum ("circular") Rydberg atoms. 

As distinct from the "linear" atoms [1, 0] , it does not 
show any sharp (threshold-like) onset of the diffusion 
with increase in the magnetic field. This is not surpris- 
ing because a distinctive feature of the linear atoms is 
singularity of their trajectories near the Coulomb center. 
Such trajectories are strongly affected by the external 
perturbations, resulting in the kink-like peculiarities (c/. 
Fig. 2 in Ref. [2]), etc. On the other hand, the electron 
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trajectories of "circular" Rydberg atoms are always far 
away from the Coulomb center and, therefore, everything 
changes smoothly. 

An important practical conclusion following from our 
consideration, e.g. for the experiments with antihydro- 
gcn, is that it might be unreasonable to increase the 
magnetic field strength above some critical value: other- 
wise, the magnetically-stimulated diffusion of the created 
atoms will override trapping the charged particles by the 
magnetic field. 

In summary, we derived an explicit analytical formula 
for the coefficient of magnetically-stimulated diffusion for 
the high-angular-momentum ("circular") Rydberg atoms 
(in contract to the previous studies, which involved nu- 
merical methods and dealt mostly with "linear" atoms). 
The case considered in our work should be well relevant 
to the Rydberg atoms formed by the recombination of 
ultracold plasmas in MOT as well as to the experiments 
on the production of antihydrogen. 

A considerable part of the present work was carried 
out during my visit to the Max Planck Institute for the 
Physics of Complex Systems (Dresden, Germany). I am 
grateful to Prof. H. Kantz for valuable discussions and 
comments. 
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